Absence of a metallic phase in charge-neutral graphene with a random gap 
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It is known that fluctuations in the electrostatic potential allow for metallic conduction (nonzero 
conductivity in the limit of an infinite system) if the carriers form a single species of massless two- 
dimensional Dirac fermions. A nonzero uniform mass M opens up an excitation gap, localizing all 
states at the Dirac point of charge neutrality. Here we investigate numerically whether fluctuations 
5M S> 7^ in the mass can have a similar efi^ect as potential fluctuations, allowing for metallic 
conduction at the Dirac point. Our negative conclusion conflrms earlier expectations, but does not 
support the recently predicted metallic phase in a random-gap model of graphene. 
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Two-dimensional Anderson localization in the Dirac 
equation shows a much richer phase diagram than in the 
Schrodinger equation.^ The discovery of graphene^ has 
provided a laboratory for the exploration of this phase 
diagram and renewed the interest in the transport prop- 
erties of Dirac fermions.^ One of the discoveries resulting 
from these recent investigationsir— was that electrostatic 
potential fluctuations V{r) induce a logarithmic growth 
of the conductivity ct oc InL with increasing system size 
L. In contrast, in the Schrodinger equation all states 
are localized by sufficiently strong potential fluctuations^ 
and the conductivity decays exponentially with L. 

Localized states appear in graphene if the carriers ac- 
quire a mass M{r), for example due to the presence of a 
sublattice symmetry breaking substrate^i^ or due to ad- 
sorption of atomic hydrogen .iiSiii Anderson localization 
due to the combination of (long-range) spatial fluctua- 
tions in M{r) and V{r) appears in the same way as in 
the quantum Hall effect (QHE)>iii^ All states are local- 
ized except on a phase boundarj*i^ of zero average mass 
M = 0, where a takes on a scale invariant value of the or- 
der of the conductance quantum Go = 4e^/ft (the factor 
of four accounts for the two-fold spin and valley degen- 
eracies in graphene). 

An altogether different phase diagram may result if 
only the mass fluctuates, at constant electrostatic poten- 
tial tuned to the charge neutrality point (Dirac point, 
at energy E = Qi). The universality class is now differ- 
ent from the QHE, because of the particle-hole symmetry 
axH*ax = —H oi the single- valley Dirac Hamiltonian 
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This symmetry is broken by an additional electrostatic 
potential, as well as by intervalley scattering. Anderson 
localization in the presence of particle- hole symmetry has 
been studied extensivelyi^^J^ in the context of supercon- 
ductivity, where the Dirac spectrum appears from the 
superconducting order parameter rather than from the 
band structure. The (numerical) models used in those 
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FIG. 1: Contour plot of a random mass with Gaussian corre- 
lator ([2]), for Kq — 10. The zero-mass contours are indicated 
in black. 



studies contain randomly distributed vortices in the or- 
der parameter, and are therefore not appropriate models 
for graphene. 

It is the purpose of this work to identify, by numeri- 
cal simulation, what is the phase diagram of the Dirac 
Hamiltonian with a random mass M(r) = M + SM{r) 
— in the absence of any other source of disorder. This 
study was motivated by recent analytical work by Ziegler 
in the context of graphene which predicted a transi- 
tion into a metallic phase upon increasing the disorder 
strength 6M at constant average mass M 7^ 0. Such 
a metal-insulator transition was known in the context 
of superconductivity)^ but it was understood that this 
requires vortex disorder 1^"— In order to resolve this con- 
troversy, we perform a numerical scaling analysis of the 
conductivity and find no metallic phase as we increase 
SM. 

We calculate the conductivity a for a two-dimensional 
strip geometry between electron reservoirs (at x = 
and X = L, see inset in Fig. [2]), with periodic bound- 



2 




FIG. 2: Average conductivity a as a, function of length L 
(for fixed W — 800 f). The average mass is set at M = 0, 
while the mass fluctuations are varied by varying Ko- The 
dashed line is at ao/Go — I/tt. The inset shows the layout of 
the disordered charge neutral strip (dotted rectangle) between 
infinitely doped electron reservoirs at a voltage difference V 
(gray rectangles). 



ary conditions in the transverse direction (at y = and 
y = W). The Fermi level is tuned to the Dirac point 
in the strip, while it lies infinitely far above the Dirac 
point in the reservoirs. For zero mass M and large as- 
pect ratio W/L the conductivity has the scale indepen- 
dent value2^i2^ ctq = Gq/tt. We generate a random mass 
with Gaussian correlator 

(W/(r.)^Af(r')) = ^^^e-l-'lV^C^ (2) 

characterized by a correlation length ^ and a dimension- 
less strength 



Ko = (v/ny / dr{SM{0)SM{r)). 



(3) 



A contour plot for a single realization of the disorder is 
shown in Fig. [TJ 

The N X N transmission matrix t through the strip is 
calculated from iJoirac by application of the numerical 
method of Ref. to a random mass rather than to a 
random scalar potential. We obtain t from the transfer 
matrix T, which relates {i/jix = L)) — T\i^{x = 0)) and 
is given by 



d_ 

dy 



-May. (4) 



Scattering from the fluctuating mass SM{r) in the slice 
{n ~ l)6x < X < nSx, of incremental length Sx = L/N^, 
is approximated by the transfer matrix 



srn = 



1 - UM„{y)a, 



V 



\5Mn{y)(Ty' 

p7l6x 



5Mn{y) = - 



{n — l)5x 



dxSM{r). 
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FIG. 3: Same as Fig. [21 but now for a nonzero average mass 
Af = 5 • IQ-'^h/v^ (solid curves, W = 800 C) and M = 5 ■ 
10"^ h/v^ (dashed curves, W — 400^). The lower panel shows 
the same data on a logarithmic horizontal scale, rescaled by 
6oc=C//(ifo,M). 



The approximation ([5]) becomes exact in the limit 
oo. Moreover, for any Nl it satisfies the requirements of 
particle-hole symmetry {axT*<Jx = T) as well as current 
conservation {axT^Ux = T~^). 

We thus obtain the conductance G = GqTy tt^ and the 
conductivity a = G x L/W. The number of transverse 
modes N and longitudinal slices are truncated at a 
finite value, which is increased until a sample specific 
convergence is reached. For the data presented, this is 
typically achieved when N = 400—800 and Nl = 300— 
600, the larger values needed for larger values of Kq. The 
sample width W = 400^ — 800^ is chosen large enough 
that the conductivity is independent of the ratio W/L. 
(Typically, W/L > 3 — 5, with the larger values needed 
for smaller values of M.) Averages over a large number 
of disorder configurations (typically 1000) produce the 
results plotted in Figs. [2] and [S] 

For M = (Fig. [2]) the conductivity stays close to 
the scale invariant value CTq (dashed line), no matter how 
large the disorder strength, while for nonzero M (Fig. 
[3]) the conductivity decays with increasing L. For suf- 
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ficiently large L/^ we expect single-parameter scaling, 
meaning that the data for different Kq and M should all 
fall on a single curve upon rescaling L f{KQ,M)L. 
(This amounts to a horizontal displacement of data sets 
on a logarithmic horizontal scale.) The length ^loc — 
^/ f{Ko, M) can then be identified with the localization 
length (up to a multiplicative constant). As one can see 
in the lower panel of Fig. [31 the data sets collapse rea- 
sonably well onto a single curve upon rescaling. (The 
remaining deviations may well be due to finite-size ef- 
fects.) 

For weak disorder {Kq < 1) our results are similar 
to earlier work on the superconducting random mass 
model.— That model however shows a metal-insulator 
transition at values of Kq — Kc of order unity 
(weakly dependent on M), such that for larger disorder 
the conductivity increases logarithmically with system 
sizei^ii^ 

a = tro ln(L/0, for K„>K,^l. (6) 

As argued by Read, Green, and Ludwig2£i2^ and by Boc- 
quet, Serban, and Zirnbauer;-^> metallic conduction in 
a random mass landscape requires resonant transmission 
through contours of zero mass (the black contours in Fig. 
[1]). These contours support a bound state at zero energy, 
if and only if they enclose an odd number of vortices. 
Without vortices, the phase shift accumulated upon cir- 
culating once along a zero-mass contour equals tt — so 
there can be no bound state and hence no resonant trans- 
mission. (The TT phase shift is the Berry phase of the 
rotating pseudospin cr in i?Dirac, without any dynamical 
phase shift because the energy is zero.) Our numerical 
finding that there is no metallic conduction in the ran- 
dom mass landscape without vortex disorder is therefore 



consistent with these analytical considerations. 

From the more recent analytical work by Ziegler— we 
would expect a transition into a phase with a scale in- 
variant conductivity 

a, = ao[l-{MlM,f], (7) 

when Mc — (?i/w^) exp(— 7r/ifo) becomes larger than M 
with increasing disorder strength Kq . The corresponding 
critical disorder strength = 7r/ln \v£_/hM\ « 0.6 — 1.0 
for the values of M in Fig. [3l The numerical findings 
of Fig. [21 with a decaying conductivity for Kq > lOi^c, 
do not support this prediction of a nonzero M^. Note 
that the numerical data of Fig. [21 with a scale invariant 
conductivity Uc = (To for M = 0, does agree with Eq. ((Tj) 
— it is the M > data that is in disagreeement. 

In conclusion, we have presented numerical calcula- 
tions that demonstrate the absence of metallic conduc- 
tion for the Dirac Hamiltonian ([1]), in a random mass 
landscape with nonzero average and dimensionless vari- 
ance Kq ^ 1. The decay of the conductivity with sys- 
tem size L is slower for larger disorder strengths, but 
no metal-insulator transition is observed. A transition 
into a metallic phase (with ct oc In L) has been attributed 
to vortex disorder 32r— Our numerical results are consis- 
tent with this attribution, since our model contains no 
vortices and has no metallic phase even if Kq 3> 1. 
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thanks the Dahlem Center at FU Berlin for hospitality 
during the inital phase of this project. 
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